Polarization states of gravitational waves with a 
massive graviton 



Wayne L.S. de Paula, Oswaldo D. Miranda and Rubens M. 
Marinho 

Institute Tccnologico de Aeronautica, Departamento de Fi'sica, Praga Marechal 
Eduardo Gomes 50, Sao Jose dos Campos, 12228-900 SP-Brazil 



Abstract. Using the Newman— Penrose formalism, we obtain the expHcit 
expressions for the polarization modes of weak, plane gravitational waves with 
a massive graviton. Our analysis is restricted for a specific bimetric theory whose 
term of mass, for the graviton, appears as an effective extra contribution to the 
stress-energy tensor. We obtain for such kind of theory that the extra states 
of polarization have amplitude several orders of magnitude smaller than the 
polarizations purely general relativity (GR), /i+ and /ix, in the VIRGO— LIGO 
frequency band. This result appears using the best limit to the graviton mass 
inferred from solar system observations and if we consider that all the components 
of the metric perturbation have the same amplitude h. However, if we consider 
low frequency gravitational waves (e.g., /gw~10~^Hz), the extra polarization 
states produce similar Newman— Penrose amplitudes that the polarization states 
purely GR. This particular characteristic of the bimetric theory studied here could 
be used, for example, to directly impose limits on the mass of the graviton from 
future experiments that study the cosmic microwave background (CMB). 
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1. Introduction 

General Relativity (GR) assumes that gravitational forces are propagated by a 
massless graviton. However, the present experimental limits on the mass of the 
graviton are only based on the behaviour of static gravitational fields as, for example, 
the Newtonian planetary motion in the solar system. If gravity is described by a 
massive graviton, the Newtonian potential would have Yukawa modifications of the 
form 

CM 

yW = — cxp(-r/Ag), (1) 

where M and Ag are the mass of the source and the Compton wavelength of the 
graviton, respectively. 

The best bound on the graviton mass from planetary motion surveys is obtained 
by using Kepler's third law to compare the orbits of earth and mars, yielding 
TTig < 4.4 X lO^^^eV 0. Another bound on the graviton mass can be established by 
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considering the motions of galaxies in clusters of galaxies, yielding < 2 x lO^^^eV. 
This second bound is less robust than solar system estimates due to the uncertainty 
about the matter distribution of the universe on large scales ^ 12 E] . 

A graviton with nonzero mass would produce several effects in the dynamical 
regime as, for example, extra degrees of polarization for the generation of gravitational 
waves and velocities of propagation dependent on the frequency of the waves. 

Based on these characteristics ^ has suggested that the mass of the graviton 
could be bounded using gravitational wave observations. As binary systems evolve, 
they will slowly spiral together due to the emission of gravitational radiation. Over the 
course of time, the frequency of the binary orbit rises, ramping up rapidly in the late 
stages of the evolution, just prior to coalescence. Laser interferometer gravitational 
wave detectors should be able to track the binary system's evolution, obtaining the 
detailed time-dependent waveform using the matched filtering techniques required for 
data analysis in these detectors. 

At least in principle, another possibility to identify the effects produced by 
massive gravitons consists in studying the excited vibrational eigenmodes of spherical 
gravitational wave detectors to identify the field content of a specific gravitational 
theory by the observed features of the waves [Ij . 

It is worth stressing that gravitational wave detectors, either interferometers or 
resonant (cryogenic bars and spheres) are members of a network that will permit the 
reduction of spurious signals and an experimental determination of the false alarm 
rate. In addition, three or more detectors ensure the complete reconstruction of a 
gravitational wave event, including the determination of its velocity of propagation 
and, the identification of the Riemann tensor signatures. 

From the theoretical point of view, massive gravitons have a number of strange 
properties produced at the level of the field equations. If the mass term has a 
specific Fierz— Pauli structure the propagator around flat space suffers the van 

Dam— Veltman— Zakharov (vDVZ) discontinuity [3011]. 

On the other hand, the effects of Fierz— Pauli mass terms were considered only 
around flat space backgrounds. Thus, the discontinuity that was found could be a 
peculiarity of the flat background and some of the difhculties could be evaded by 
considering a background with curvature jlO|. 

In fact, it was found that in constant curvature backgrounds (AdS spaces), the 
extra polarizations of the massive gravitons have a coupling m^/ H where H is the 
Hubble constant of the AdS space |11II12| . In this case, the predictions of the massive 
theory were the same as of the massless theory when mg <C H . 

Nevertheless, it is important to stress that GR has an excellent agreement in the 
prediction of the decrease of the orbital period (r) of the binary pulsar PSR 1913+16 
at least in the weak field limit. The GR predicts r = —2.4 x 10~^^ and the observed 
value is T = -(2.40243 ± 0.00005) x lO'^^ j^j. Thus, an alternative theory of gravity 
in the limit mg should obtain the same results of GR. Besides, the Newtonian 
limit has to be valid. 

In fact, theories of gravitation can be divided into two different groups: metric 
and non-metric theories |T3] . Basically, a theory of gravitation is said to be metric if 
it obeys the principle of equivalence. That is, the action of gravitation on the matter 
is due exclusively to the metric tensor. GR and Brans— Dicke theories are examples 
of metric theories of gravitation. 

Thus, one of the ways to obtain a massive graviton theory, preserving the 
equivalence principle, is to add a prior geometry. This possibility was recently 
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proposed by |14| . In his theory, it is possible to recover GR when mg — > without the 
problem of the vDVZ discontinuity because, the linearized mass term of the theory is 
not a Pauli— Fierz term. 

In the present study, we analyse the polarization states of the gravitational waves, 
in the theory developed by using the formalism developed by Newman— Penrose 
[T^ . Then, we follow the method proposed by [Ifil [T7] to determine the number of 
polarization modes of the gravitational waves. The final result consists in the explicit 
expressions for the six independent 'eletric' components of the Riemann tensor. 

In section 12 we present the field equations for the massive graviton using 
the action proposed by In section 13 we obtain the general solution for 

weak gravitational waves with a massive graviton. In section 01 we determine the 
explicit expressions for the polarization modes of the gravitational waves using the 
Newman— Penrose formalism. In section we present our conclusions and discuss 
some aspects on how it could be possible to impose limits on the mass of the graviton 
from observations of the cosmic microwave background. 



2. Field equations for massive graviton 

We wish to examine an extension of linearized GR which includes a massive graviton in 
the field equations. As exposed above, our main idea is recover GR when mg — ^ so, 
the linearized mass term is not the Pauli— Fierz term. That is an essential condition 
to have a well-behaved classical limit as the graviton mass goes to zero. According to 
| 14| a way to get this is to add the massive graviton term in the action as an isolated 
term. That is 

r C /q /p ^mass ,n\ 

IdttG c c c 

The last term on the right-hand side of (O contributes to the equations of motion 
with an effective stress tensor given by 

TA'J' — ('Q^ 

-'-mass ~ , r -'mass, 1,0; 

where 



1 



(5o"': 



(3-5or 



igorigo'Uig - 9or^ 



(do ^)pi^ 



(4) 



where (50)0/3 is the non-dynamical background metric, gafj is the physical (dynamical) 
metric and mg is the graviton mass. 

The weak field limit is obtained with g^p = igo)ai3 + ha/3 and {hapl ^ |(5o)a/3|- 
In reality, as discussed by |14|. any algebraic function of the physical metric and 
background metric with correct linearized behaviour up to second order in h would 
present the same characteristics of equation 



Then, in the weak field limit, 



T^J''''' becomes 



rpmass 



{go'r^haf3 



(go] 



fJ.V 



(5) 
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The field equations can be rearranged to produce a structure like the usual 
Einstein field equations. That is, 

(6) 

where T^" represents the usual energy-momentum tensor. 

From the last equation it is possible to verify that when m,^ 0, we recover the 
ordinary Einstein field equations. 



3. Weak gravitational waves with a massive graviton 

Considering gravitational waves far from field sources, we can work in the weak-field 
limit. In this case, we should have to choose a particular background geometry for 
the non-dynamical metric. The natural choice, in a first work, is to take (50)^11^ to 
correspond to a flat spacetime (Minkowski spacetime). In principle, this choice should 
produce a good agreement with all astrophysical observations at the level of weak 
gravitational fields. Thus, 



(7) 



where ^ 1, is the Minkowski spacetime metric and the metric signature is 
(-!,+!, -M,+l). 

Far from gravitational sources, we can take T^"^ = (vacuum solution) and we 
have from the conservation of stress-energy 



\7 T^M'^ 

4 ^ M mass' 



where G^"^ is the Einstein tensor and V denotes the covariant derivative. 
Thus, we have for the right-hand side of equation 



ft2 



1 



(9) 



where we used in the last equation the result obtained from equation 0. 

It is worth stressing that equation © is not a gauge choice but a constraint 
imposed by the conservation of energy. 

On the other hand, the Ricci tensor is to first order in h 



R.J 



d 



■pa 



d 



and the affine connection is 

dxt^'^"" ' dx'' 
Thus, equation (|SJ) can be rewritten as 
1 



9 u 9 r 



d 



□ h 



where we defined = rn^ (? jY? . 

The last equation rewritten in terms of the trace reverse of h^y produces 



0, 



(10) 



(11) 



(12) 



(13) 
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where h^j,^ = /i^^ - l/2ry^^/i. 

The general solution of equation is a linear superposition of solutions of the 

form 

hf,^ = ef,^exp{ikax"^), (14) 

where e^^ is the polarization tensor. 

From the normalization condition k^^k^ = —rn^, we can obtain the dispersion 
relation 

k = V(c^/c)2 -m2. (15) 

As a consequence of massive gravitons, the speed of propagation of a gravitational 
wave is dependent on the frequency and given by v[lj) = duj/dk which produces 

v{.) = c^l-^. (16) 

It is important to stress that we are considering in our study a gravitational wave 
travelling in the +2:-direction. Thus, the wave vector is given hy kz = k and the metric 
perturbation presented in equation (14) may be rewritten as 

^/ii/ = e^jy exp(— i t + z A: z). (17) 

In the next section we will introduce the tetrad formalism to determine the 
polarization wave modes in this massive theory. 



4. Polarization states of gravitational waves v^rith a massive graviton 

For nearly null gravitational waves in the weak field limit, it is only necessary to restrict 
our study to the form and behaviour of the Riemann tensor. It is the Riemann tensor 
that gives us the information how a gravitational wave interacts with a detector. 

To analyse the components of the Riemann tensor into independent wave modes in 
as invariant a manner as possible, we should investigate the transformation properties 
of the Riemann tensor under Lorentz transformations which leave the wave direction 
fixed. 

Then, we choose a basis vectors in which the components of the Riemann tensor 
are computed. The basis vectors form a quasiorthonormal tetrad basis. In particular, 
we follow the formalism derived by 15 and used in the work of ^16,. .17 . 

It is worth stressing that the original tetrad formalism differs from the 
Newman— Penrose formalism only in the manner of choice of the basis vectors. That is, 
instead of an orthonormal basis, the choice is made of a complex null-basis (fc, I, m, m) 
where k and I are two real null-vectors and m and fh are a pair of complex conjugate 
null-vectors. 

These vectors satisfy the orthogonality relations 

k ■ I — I, m ■ fh = —1, k ■ m — k ■ fh — I ■ m — I ■ fh — 0. (18) 

In particular, k and I are tangent to the propagation directions of the two plane 
waves, a wave propagating in the +z direction and a wave propagating in the —z 
direction, respectively. 
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Following |16l I17| and taking into account the normalization conditions in 
equation (|18|) . the vectors {k, I, m, m) can be written as 

fc = --^(1,0,0,1), (19) 
/ = --^(1,0,0,-1), (20) 
m= --^(0,l,i,0), (21) 

m = -^(0,l,-i,0). (22) 

Although this basis vector form a null-basis, it is possible to expand the wave 
vector k for a gravitational wave not exactly null ^3], like that obtained from the 
bimetric theory studied here. 

In the Newman— Penrose formalism, the Riemann tensor is split into irreducible 
parts: the Weyl tensor, the traceless Ricci tensor and the Ricci scalar named, 
respectively, tetrad components <i>, and A. 

The ten independent components of the Weyl tensor are represented by the five 
complex scalars, 

I'd = -Cpgrsk" m« k'' , (23) 
*! = -CpqrskP l"^ ¥ , (24) 
1'2 = -CpqrskP TO« fh'' P , (25) 

= -CpgrskP Z« m'' P , (26) 
-Cpqrsl^m'^l''m' . (27) 
It is also possible to define the following scalars representing the Ricci tensor 



and 



$00 = ^ 2 ^'^^ ' 


(28) 


$01 — $10 = ^ 2 ^^"^ ' 


(29) 


$11 = — ^ Rkl + Rmfh , 


(30) 


$12 = $21 - -^Rlm , 


(31) 


$22 = —-^Rll , 


(32) 


$02 = $20 = ~ 2 ' 


(33) 


R {Rkl — Rmfh) 

^ 24 ~ 12 


(34) 



In GR where v = c only the component ^'4 is not null. However, in the bimetric 
theory studied here, we observe that for gravitational waves with v{uj) ^ c, the tetrad 
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components are *o - ©(e^i?), ^'i ~ 0(ei?), $00 ^ 0{e'^ R), $01 - *io - $02 - 
$20 — 0(ei?), $11 = 3/2 *I'i2 — where e is related to the difference in speed 
between light and the propagating gravitational wave. That is, e = (c/w(u;))2 — 1. 

Thus, to describe the six independent components of Riemann tensor we shall 
choose the set \l/2, ^'3, 5*4 and $22- 

It is important to stress that for low frequency gravitational waves where e > 1, 
the componentes \l/o, vE'i, <&oo, 'I'oi, ^'107 $02, *i>20 can be written in function of 
4*2, ^3, ^'4 and <I>22- Thus, these six independent components are able to describe 
completely the polarization modes of a gravitational wave. 

Then, the explicit expressions for these six components are 

= 12 iu^TW) ' 

1 /ii3(w^ fc - cjfc^ - fc''') + i/i23(w3 + Cc)2 fc - tjfc2 - , ^, 
*3 = o , (36) 



*4 = \[hm{-Lo'^ -2u)^k + 2ujk^ + k^) 
8 



- /l22(2 + fc + 4^2 ^2 ^ fc3 _^ 2 

+ hzzi-uj^ -2Lo^k + 2ujk^ + k^) 

+ ihi2{2Lo^ + 'iuj^k + 'iuj'^k^ + Aujk^ + 2 fc'*)]/(w2 + k^), (37) 



1 (fepo + /i33)(-^' ~ 2 u;^ fc + 2 fc3 + fc4) 

*2^-8 ^:;^tp^ ' 

where we use in the above set of equations c = 1. 

In figure 1 we show the displacement that each mode induces on a sphere of test 
particles. The wave propagates in the +z direction and has time dependence cos oj t. 
The solid line is an instantaneous at cijt — and the broken line is one at tot = n. 

We can solve the set of equations , iHI^ , llT7|) , for determined gravitational 
wave frequencies, in order to identify the contribution of the extra polarization modes 
for a signal received by a gravitational wave detector. 

In particular, we present the results for three different frequencies: /gw ~ 100 Hz 
that corresponds approximately to the frequency of maximum sensibility for the LIGO 
interferometer; /gw = 10"'^ Hz that corresponds to the maximum sensibility for the 
future laser interferometer space antenna (LISA) and /gw — 10^^ Hz. 

Thus, we have rUg — 0.44 x 10^2icV/c2 that corresponds to the best limit from 
solar system observations jl8| . 

a) /gw = 100 Hz: 

*2-2.1x 10-35(/^33_/^oo), (39) 

*3- 1.2X 10-3^(/ll3+i/l23), (40) 

*4 ^ 1.2 X IQ-^^i-hno - 1133) + 4.4 x 10-iS(-/z22 +1/112), (41) 

$22^1.2X 10-34(_/joo_/j33). (42) 

b) /GW = 10-3Hz: 

*2^2.1x 10-35(/i33_/joo), (43) 
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Figure 1. The six polarization modes of weak gravitational waves permitted in 
any theory of gravity. Shown is the displacement that each mode induces on a 
sphere of test particles. The wave propagates out of the plane in (a), (b), (c), and 
it propagates in the plane in (d), (e) and (f). The displacement induced on the 
sphere of test particles corresponds to the following Newman— Penrose quantities: 
Re *4 (a), Im *4 (b), #22 (c), ^2 (d), Re *3 (e), Im *3 (f)- 



*3 - 1.2 X 10-34(/j^3+ 1/^23), (44) 

*4 ^ 1.2 X 10-34(-/ioo - /133) + 4.4 X 10-26(-/i22 + i/ii2), (45) 

$22 ^ 1.2 X lQ-^\-hoo - /133). (46) 
c) = 1.1 X 10-^ Hz: 

*2 ^ -2.1 X lO-^^/ioo + 3.9 X 10-^^/133, (47) 

*3 ^ 7.6 X 10-35(/ii3 + i/123), (48) 

*4 ^ 8.9 X 10-3^(-/ioo - /133) + 2.0 X 10-3^(-/i22 + i/112), (49) 

$22 ^ 8.9 X 10-3^(-/ioo - /133). (50) 
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5. Final remarks 

Using tlie bimetric tlieory proposed by jl4| . we obtain the explicit expressions for the 
polarization modes of gravitational waves considering a nonzero mass for the graviton. 
In particular, we analyse what happens to the six tetrad independent components for 
three different frequencies, namely, /gw^ = 100 Hz approximately the frequency of 
maximum sensibility for the LIGO interferometer; /gw = 10^^ Hz the frequency of 
maximum sensibility for the future laser interferometer space antenna (LISA) and 
/gw^ lO-'Hz. 

From the above results, we can see that for /gw — 100 Hz the dominant tetrad 
component is 4*4, with the part purely GR several orders of magnitude bigger than 
the extra polarization terms. For example, ^'4 can be split into 'I'4 ~ (^'4)gr + (^4)m, 
where the polarizations purely GR are given by (5'4)gr — 4.4 x 10^^^(— /122 + i/112) 
while, the term due to the nonzero mass for the graviton is ('I'4)m — 1.2 x 10^'^'*(— /iqo — 

/l33)- 

If we consider that all metric perturbations /i^^ have a similar value h then we 
obtain |(^4)gr| — 10^* |(^4)m|- Thus, the term (^'4)m has behaviour similar to a very 
small perturbation when compared to the GR term in ^4. A similar conclusion can be 
obtained from the comparison of the amplitudes of gravitational waves in equations 
(jSni, lEOI), (EH and (03). We can see that the modes ^'2, *3 and <i>22 maintain the 
same very small perturbation behaviour as the component (5'4)m- 

Then, at least for rUg — 0.44 x lO^^^eV/c^ and taking into account the weak field 
limit studied in the present paper, it should not be possible to identify the signature 
of the extra polarization modes from an astrophysical source detected, for example, 
by VIRGO and LIGO interferometers. Certainly, as mentioned above, this conclusion 
is based on the fact that all metric perturbations h^^ have a similar value h. In this 
case, the extra polarization modes are very small perturbations of (^4)gr- 

For the case /gw = 10^"^ Hz, we can see that the extra polarization modes 
have the same amplitudes as in the frequency /gw = 100 Hz. However, the 
weight of the term purely GR in ^"4 decreases when compared to the massive term 
(^4)m- For example, from equation (|45|l we obtain that |(^'4)gr| — 10^|(5'4)m| for 
/gw = 10-3Hz. 

As a consequence, we could think that, in principle, a gravitational wave signal in 
the LISA frequency range could give us more information on the polarization modes 
\E'2, ^'3 and $22 than a gravitational wave signal in the VIRGO-LIGO frequency 
range. But unfortunately, same as for /gw = 10^^ Hz the extra polarization modes 
correspond to very small perturbations of the GR polarizations. 

On the other hand, a very interesting result appears for the frequency /gw = 
1.1 x 10" ''Hz. In this case, we can see from equations (|T7|) . (gHJ and (|55|l 

that all tetrad components have similar amplitude. As a consequence, we have 
|('I'4)gr| — |(^4)m|- In particular, all the polarization modes produce similar 
'excitation' in the frequency 1.1 x 10~^Hz. 

This kind of result is maintained if we change the value of the graviton mass. 
However, if we reduce the value of the graviton mass, the main result is to shift the 
frequency /c where all the polarization modes have similar amplitude. That is, lower 
the graviton mass, lower the frequency /c is. 

In figure 2 we can see the relation between the mass of the graviton and the 
frequency /c. For example, if the mass of the graviton is mg — 0.44 x 10~^^eV/c^, 
then we have /c ~ 10"^ Hz. On the other hand, if the mass of the graviton is 
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'I ' '"I ' '"I ' '"I ' '"I ' '"I ' '"I ' '"I ' '"I ' '"I ' '"I ' '"I ' '"I ' '"I ' '"I 




-iQ-36 r 

^o^^ 10" 1x10" 1x10" 1x10' 

4 (Hz) 



Figure 2. The relation between tlie mass of tlie graviton and tlie frequency where 
all the polarization modes presented in figure 1 have the same amplitude. 



mg = 0.44 X 10-29eV/c2 then /c ~ lO^i^Hz. 

This resuh could be used to obtain a new hmit on the graviton mass based on 
the analysis of the maps produced by experiments that study the anisotropy of the 
cosmic microwave background (CMB). 

In particular, the generation of a stochastic background of primordial 
gravitational waves is a fundamental prediction of inflationary models for the early 
universe. Its amplitude is determined by the energy scale of inflation, which can widely 
vary between different inflationary models. 

Gravitational wave detectors, however, are quite unlikely to have enough 
sensitivity to detect such a primordial signal, owing both to its smallness and to 
its extremely low characteristic frequencies. The existence of ultra-low-frequency 
gravitational radiation, however, can be indirectly probed thanks to the temperature 
anisotropy and polarization it induces on the CMB radiation. 

In particular, the curl component, called B-mode, of the CMB polarization 
provides a unique opportunity to disentangle the effect of tensor (gravitational- wave) 
from scalar perturbations, as this is only excited by either tensor or vector modes 

From this point of view, future satellite missions, such as Planck, which will 
have enough sensitivity to either detect or constrain the B-mode CMB polarization 
predicted by the simplest inflationary models, might represent the first space-based 
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gravitational- wave detector . 

Thus, future CMB missions could present an alternative way to impose a new 
upper limit on the mass of the graviton and to constrain the number of polarization 
modes of the gravitational waves. 
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